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Abstract 

We investigate whether the anomalous tensor operators with the Lorentz structure 
+ 75) ® <7 M1/ (1 + 75), which could provide a simple resolution to the polarization 
anomaly observed in B — > (fiK* decays, could also provide a coherent resolution to the 
large B(B — > 17K*) and survive bounds from B — ► rjK decays. Parameter spaces satisfying 
all these experimental data are obtained, and found to be dominated by the color-octet 
tensor operator contribution. Constraints for the equivalent solution with (1 +75) 0(1+75) 
operators are also derived and found to be dominated by the color-singlet one. With the 
constrained parameter spaces, we finally give predictions for B s — ► (jxp decay, which could 
be tested at the Fermilab Tevatron and the LHC-b experiments. 

PACS Numbers: 13.25.Hw, 12.38.Bx, 12.15.Mm. 
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1 Introduction 



Looking for signals of physics beyond the Standard Model (SM) is one of the most important 
missions of high energy physics. Complementary to direct searches for new physics (NP) parti- 
cles in the high energy colliders, the study of B physics is of great importance for probing indi- 
rect signals of NP. In this respect, the B factories at SLAC and KEK are doing a commendable 
job by providing us with a huge amount of data on various 5-meson decays, which are mostly 
in perfect agreement with the SM predictions. However, there still exist some unexplained 
puzzles, such as the unmatched CP asymmetries in B — > ttK decays [HE], [3], the abnormally 
large branching ratios of B — > rj'K and B — > r\K* decays [UI21I11IH1E1EIE1IH1 > an ^ the l &r g e 
transverse polarization fractions in B -> <f)K* decays [H [21 [HI H21 LI31 [HI LIS1 [HI UHl LTHl M- 
Confronted with these anomalies, we are forced not only to consider more precise QCD effects, 
but also to speculate on the existence of possible NP scenarios beyond the SM. 

It is well-known that the flavor-changing neutral current (FCNC) processes arise only from 
loop effects within the SM, and are therefore very sensitive to various NP effects. Since the 
puzzling .B-meson decay processes mentioned above are all related to the FCNC b — > s transi- 
tions, these decay channels could be used as effective probes of possible NP scenarios. So, if one 
kind of NP could resolve one of these puzzles, it is necessary to investigate whether the same 
scenario can also provide a simultaneous resolution to the others. Considering the fact that 
current theoretical estimations of BiB — > rj'K) still suffer from large uncertainties j9j [21], and 
the NP scenario with anomalous tensor operators considered in this paper do not contribute 
to the B — > (pK and B — > irK decays in the naive factorization (NF) approximation, we shall 
only focus on the B — > r)K^ and B — > <pK* decays. 

The recent experimental data on the longitudinal polarization fraction /x in B° — > <pK*° 
decay is given as 

0.52 ± 0.05 ± 0.02 BABAR [UJ, 
h=\ 0.45 ±0.05 ±0.02 Belle [12], (1) 

0.57 ±0.10 ±0.05 CDF P!. 
On the other hand, since the two final-state light vector mesons and K*° in this decay 
mode are flying out fleetly in the rest frame of B meson, and the structure of the charged 
weak interaction current of the SM is left-handed, as well as the fact that high-energy QCD 



2 



interactions conserve helicity, any spin flip of a fast flying quark will be suppressed by one 
power of 1/rrib, with m& the b quark mass. It is therefore expected that, within the SM, both 
of the final-state hadrons in this decay mode are mainly longitudinally polarized, with 

h ~ 1 - 0(l/m 2 b ), (2) 

while the transverse parts are suppressed by powers of m^^x* / m B- Obviously, the experimental 
data Eq. ([T]) deviates significantly from the SM prediction Eq. ([2]), and this polarization anomaly 
has attracted much interest in searching for possible theoretical explanations both within the 
SM and in various NP models [H EH E3 EH [20] . For example, the authors in Refs. pH E3 
[TUl |2"U] have studied this anomaly and found that the four-quark tensor operators of the form 
5(7^(1 +75)6 1+75)5 could offer a simple resolution to the observed polarization anomaly 
within some possible parameter spaces. 

Since the B — > r/K^ decays, in analogy with the B — > <pK* decays, also involve the b — > sss 
transition, it is necessary to investigate the effects of these new types of four-quark tensor 
operators on the latter. In particular, it is very interesting to see whether these new four-quark 
tensor operators with the same parameter spaces could also simultaneously account for the 
measured B(B — > t]K*) [21 SI El [0], which are much larger than the theoretical predictions within 
the SM [HI [9] , and survive bounds from B —> rjK decays. Motivated by these speculations, in 
this paper, we shall investigate the effects of the following two types of tensor operators on 
these decay modes (with i, j the color indices) 

T i = s<v(l + 75 )6 <g> sa^il + 75 )s, T8 = 5^(1 + 75)^- <g> 8^(1 + 75)^, (3) 

and try to find out the allowed parameter spaces characterized by the strengths and phases 
of these new tensor operators that satisfy all the experimental constraints from these decays. 
Moreover, since the (pseudo-) scalar operators 

O s+P = s(l + 75)6 ® s(l + 7 5 )s, O54-P = + T5)6j ® ^-(1 + 7 5 )si , (4) 

can be expressed, through the Fierz transformations, as linear combinations of the new tensor 
operators Eq. ([3]), constraints on these two operators can be derived easily from those on the 
latter. 
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To further test such a particular NP scenario with anomalous tensor operators, we also give 
predictions for the branching ratio and the longitudinal polarization fraction of B s — > <fxp decay, 
which is also involved the same quark level b —* sss transition. All these results could be tested 
at the Fermilab Tevatron and the LHC-b experiments. 

The paper is organized as follows. Sec. 2 is devoted to the theoretical framework. After a 
brief entertainment of the QCD factorization formalism (QCDF) [22], we discuss the anomalous 
tensor operator contributions to the B — > r/K^ and B — > <pK* decays. In Sec. 3, our numerical 
analysis and discussions are presented. Sec. 4 contains our conclusions. Appendix A recapitu- 
lates the amplitudes for the five decay modes within the SM [15j E3]- All the theoretical input 
parameters relevant to our analysis are summarized in Appendix B. 



2 Theoretical Framework 

2.1 The SM results within the QCDF framework 

In the SM, the effective Hamiltonian responsible for b — > s transitions is given as [21] 



c r 10 
v ub v: s (CiO? + c 2 o u 2 ) + VcbV* (do? + c 2 o c 2 ) - v tb v* ( £ qo, 



i=3 



-h.c, (5) 



where V q bV* s (q = u,c and t) are products of the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
elements [25] , the Wilson coefficients, and Oi the relevant four-quark operators whose explicit 
forms could be found, for example, in Ref. 



In addition to 7i c ff, we must employ a factorisation formalism of hadronic dynamics to 
study the B~ — * i]K~ , B° riK° , B~ — > r\K*~ , ~B° r]K*° , and B° <fK*° decays. To this 
end, we take the framework of QCDF [22]. The factorization formula allows us to calculate 
systemically the hadronic matrix element of operators in the effective Hamiltonian Eq. (0) 

(M 1 M 2 \O i \B) = V Ff^ Ah I dxT^x^M^x) + (M, <- M 2 ) 

■*1 rl rl 

dU dx dyT^{i,x,y)^ B {^ Ml {x)^ M2 {y), (6) 
) Jo Jo 
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where is the B — > M transition form factor, Tj- and T? 1 are the perturbatively calculable 

hard kernels, and $x( x ) (X = B, Mi j2 ) are the universal nonperturbative light-cone distribution 
amplitudes (LCDAs) of the corresponding mesons. 

In the recent years, the QCDF formalism has been employed extensively to study non- 
leptonic B decays. For example, all the decay modes considered here have been studied com- 
prehensively within the SM in Refs. [151 E3] . We recapitulate the amplitudes for B — > r/K^ 
and B — > <pK* decays in Appendix A. 

It is noted that, along with its many novel progresses in non-leptonic B decays, the frame- 
work contains estimates of some power corrections which can not be computed rigorously. These 
contributions may be numerically important for realistic 5-meson decays, especially for some 
penguin-dominated decay modes [T5"| [T7] [221 [2E]- In fact, there are no reliable methods avail- 
able at present to calculate such contributions. To give conservative theoretical predictions, at 
least, one should leave the associated parameters varying in reasonable regions to show their 
possible effects. In this work, following closely the treatment made in Refs. [221 [21], we will 
parameterize the end-point divergences associated with these power corrections as 

In the following numerical calculations, we take the parameter A^ and the phase <pA varying in 
the range 0.2 ~ 0.8 GeV and —45° ~ 45°, respectively. 

In our calculation, we have neglected possible intrinsic charm content and anomalous gluon 
couplings related to the meson 77, both of which have been shown to have only marginal effects on 
the four B — > rjK^*' decays [El [9]. As for the 77-77' mixing effects, we shall adopt the Feldmann- 
Kroll-Stech (FKS) scheme [21] as implemented in Ref. [9]. A recent study and comparison of 
different 77-77' mixing schemes has been given in Ref. [29|. 

In the amplitude for B° — >• <fiK* decay in Eq. ( 1391) . a new power-enhanced electromagnetic 
penguin contribution to the negative-helicity electroweak penguin coefficient af'sw as ^ TS ^ 
noted by Beneke, Rohrer and Yang [TS], has also been taken into account in our calculation. 
In a recent comprehensive study of B — > VV decays [TB] , it is found that the small Jl could be 
accommodated within the SM, however, with very large theoretical uncertainties. 
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2.2 Anomalous tensor operators and their contributions to the de- 
cay amplitudes 

Since the SM may have difficulties in explaining the large BiB — > rjK*) and the measured 
polarization observables in B — »• <pK* decays, we shall discuss possible NP resolutions to these 
observed discrepancies. Specifically, we shall investigate whether these discrepancies could be 
resolved by introducing two anomalous four-quark tensor operators defined by Eq. ([3]). 
We write the NP effective Hamiltonian as 



H eff ~ ^ \ V ts\ e 



+ h.c, (8) 



with the tensor operators Ot\ and Ots defined by Eq. ([3]). The coefficient Cti(ts) describes the 
relative interaction strength of the tensor operator 0ti(T8)> an d is the NP weak phase. In 
principle, such four-quark tensor operators could be produced in various NP scenarios, e.g., in 
the Minimal Supersymmetric Standard Model (MSSM) [30,131]. Interestingly, the recent study 
of radiative pion decay ir + — > e + i"~f at PIBETA detector [32] has found deviations from the 
SM predictions in the high-£' 7 -low-i? e + kinematic region, which may indicate the existence of 
anomalous tensor quark-lepton interactions [331 EH [35] . 

At first, we present the NP contributions to the decay amplitudes of B — > r/K^ and 
B — > cf)K* decays due to these new tensor operators. Since their coefficients are unknown 
parameters, for simplicity, we shall only consider the leading contributions of these tensor 
operators. 

n 



B(B-) 




(j k { *\k^-) 



d(u) 

Figure 1: Feynman diagram contributing to the decay amplitudes of B — > r/K^ decays due to 
the anomalous tensor operators. Another type of insertion has no contribution. 

As for the four B — > 7]K^ decays, the relevant Feynman diagram due to the tensor op- 
erators Oti,t& is shown in Fig. [TJ It is easy to realize that the amplitude corresponding to 
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another type of insertion would vanish in the leading order approximation. Instead of using the 
Fierz transformation, an easy way to calculate the amplitude in Fig. [1] is to use the light-cone 
projection operator of the meson rj in momentum space [23J 



75 ® v (x) - /^ s 75 $ P (x) 



(9) 

a/3 



(V(q)\sw\0) = -if°q», Ns = I£r* = (11) 



where q, § v , and $ p are the momentum, leading-twist, and twist-3 LCDAs of the meson rj, 
respectively, and k% denote the momenta of the quark and anti-quark in the meson t], and 
are given by 

&i = xq» + + <f , fc£ = xq» -h* + (10) 

2xq ■ q 2xq ■ q 

with x = 1 — x, and x the momentum fraction carried by the constituent quark. The decay 
constant /* and the factor /i^ in Eq. (jHj) are defined, respectively, by [9] 

K 

2 2f[m s 

where we have used \r)) = cos (j)\r) q ) — sin 0|?7 S ), with \r) q ) = (\uu) + \dd)) /y/2 and \r) s ) = \ss) 

After some simple calculations, the NP contributions to the decay amplitudes of the four 
B — > t)K(*> decays due to Ti-^/f i n Eq- (E]) can be written as 

^ P ^- = i ^§ 1^1 e*3^ (m% u - m\-) F^H) K > ( 12 ) 

^SL^o = t ^ |V£.| e*3»rT* (m| d - m^ ) itf^K) /* , (13) 

AT^ v k*- = -iV2G F \V ts \e iS T3g T rl°m K *-(e*-p)A^ K *(m 2 TI )fs, (14) 

^_, = -iV2G F \V u \e^3g T r%m K *(el.p)A*^(ml)f;, (15) 

where g-r = CVg + Cti/N c , and the factor 3 is due to contractions of the involved 7 matrices. It 
is interesting to note that the above four decay amplitudes are all proportional to the "chirally- 
enhanced" factor r^ s = , which has been found to be very important for charmless 

hadronic B decays [23] . 

We now present the NP contribution to the decay amplitude of 5° — > <fK*° decay. Based 
on the observation that they contribute only to the transverse polarization amplitudes but not 
to the longitudinal one [TTI [T8l [19], these anomalous four-quark tensor operators have been 
proposed to resolve the polarization anomaly observed in B — > <j)K* decays. The relevant 



s i(i)\ / Si 




bj T 1(S) S i(j) 

R° 1 I Pi 7?*° 

B U . J 7 " 

(a) 

Figure 2: Feynman diagrams contributing to the decay amplitude of B — ► ^.K" decay due to 
the anomalous tensor operators. 

Feynman diagrams are shown in Fig. For Fig. [2] (a), we shall use the following matrix 
elements [HI E3 EE] 

(0(g, £l )l<^|O) = -#(e£ft, - e£fc) , (16) 
{K*(j>',e 2 )\sa» v q»b\B(j>)) = e^e* W 2T^ K * (q 2 ) , 
(K*(p\e 2 )\rsa^ l5 b\B(p)) = (~i)T 2 B ^ (q 2 ) {e^ml - m 2 K .) - (e* 2 -p) (p + p') f 

+ H)T 3 B - K \q 2 )(el-p)U- f (P + P')A- (17) 

For Fig. [2](b), we shall use the light-cone projector operator of the transversely polarized vector 
meson [T5], IHB] 

• fT 

Mf = -^-^ ± #*±(x) + ..., (18) 

where is the leading twist LCDA of the meson <p, and the ellipsis denotes additional 

parts that have no contributions in our case. It is easy to find that neither Fig. [2(a) nor (b) 
contributes to the longitudinal polarization amplitude, and the final decay amplitude can be 
written as 

^ = ^|^|e^^(-4z/J)|^ pCT ^4V^2T^*(mJ) 
+T 2 M *K) [(e* 1± ■ el) (ml - m 2 K «) - 2(e* 1± ■ p)(s* 2 ■ p) 
-22f^(mJ) J (sl ± -pM-p) , (19) 
with g' T = (1 + ^r) Ctx + (| + jf-) C*ts- in the helicity basis, the new decay amplitude Eq. (|TU1) 
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can be further decomposed into 

<o P = 0, (20) 



o-NP _ Gf , t/ I iS T / ( A - f T\ 



[m\ d - m 2 K ,o) T^ K \ml) T 2m Bd p c if-** (mf 



(21) 



where p c is the center-of-mass momentum of final mesons in B° rest frame. Compared with 
the SM predictions Eqs. (I40p and (|4ip . the new transverse polarization amplitudes H±± are 
enhanced by a factor of m^/m^, while the longitudinal part remains unchanged. It is therefore 
expected that these new tensor operators might provide a possible resolution to the polarization 
anomaly observed in B — > cf)K* decays. 

2.3 The branching ratios and polarization fractions 

From the above discussions, the total decay amplitudes are then given as 

^ = ^ SM + ^ NP , (22) 

where ^4 SM denotes the SM results obtained using the QCDF, and „4 NP the contributions of 
the particular NP scenario with anomalous tensor operators in Eq. (jSJ). The corresponding 
branching ratios are 

B(B°> + -> V K^°' + ) = I^L\A(B°< + -> v K^ ' + )\ 2 } (23) 

8irm B 

B&^jF ) = ^(\Hoo\ 2 + \H ++ \ 2 + \H^\ 2 ), (24) 

where tb is the life time of B meson. 

In the transversity basis [39], the decay amplitude for any B — > VV decay can also be 
decomposed into another three quantities Aq, Ah, and A±, which are related to the helicity 
amplitudes H 00 , H ++ , and H through 

A = H 00 , A|, = g±±-^g= A ± = - H++ ^ H - . (25) 

In terms of these quantities, we can express the longitudinal polarization fraction as 



h 



I A) 



2 







l lll 2 + 


2 


2 


#00 


l 2 + l 


H ++ 


2 +\ 





(26) 



In addition, the relative phases between these helicity amplitudes <p»± = Arg(A^±/A ) + 
7r (the definition of these observables is compatible with that used by the BABAR and Belle 
collaborations pTTl [T2] ) . are potentially very useful for constraining the parameter spaces of NP 
scenario, however, depend on whether the strong phases of the helicity amplitudes could be 
calculated reliably. 

3 Numerical analysis and discussions 

With the theoretical formulas and the input parameters summarized in Appendix B for the 
decay modes of our concerns, we now go to our numerical analysis and discussions. 

As shown by the NP decay amplitudes Eqs. (jT2T) — (TT5T) and (l20l) - (l2TT) . the allowed regions for 
the parameters and St can be obtained from the four measured B(B°' + — > r]K^°' + ), and 
the ones for g' T from B(B° — > <f)K*°) and J'l, respectively. Generally, we have five branching 
ratios and one polarization fraction, but only three free parameters: one NP weak phase St and 
two effective coefficients gx and g' T (or equivalently Ct\ and Cts)- So, it is easy to guess that 
these decays could severely constrain or rule out the NP scenario with two anomalous tensor 
operators Eq. (E]). 

To make the guess clear, we shall first find out the allowed regions for the parameters gx 
and St from the four B°' + — > r]K^°' + decays. Then, using the allowed regions for the weak 
phase St, we try to put constraint on the parameter g' T from B(B° — > <pK*°) and fi- Finally, 
we can obtain the allowed parameter spaces, if there are, for Cti, Cts, and St that satisfy all 
the experimental data on the decay modes of our concerns. To further test the particular NP 
scenario Eq. (|BD, we also present our theoretical predictions for B s — > cfxf) decay. Theoretical 
estimations of the annihilation contributions at present suffer from very large uncertainties, 
which of course will dilute the requirement of NP very much. To show the dilution, we will 
give our numerical results for two cases, i.e., with and without annihilations for comparison. 

3.1 Constraints on the NP parameters from B° ,+ — > r]K^ 0,+ decays 

Since the tensor operator contributions to the decay amplitudes of the four B°' + —> r]K^°' + 
decays are all characterized by the parameters gT and St, possible regions for these two NP 
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parameters can be obtained from the measured branching ratios of B°' + — > r]K^°' + decays. 
Our main results are shown in Tables [TH2] and Figs. [3HU The experimental data listed in 
Table [Hare taken from the Heavy Flavor Averaging Group (HFAG) [2j. The SM predictions for 
the branching ratios of these four decays are presented in the third column of Table [TJ where 
the theoretical uncertainties are obtained by varying the input parameters within the regions 
specified in Eq. (J3) and Appendix B. For each decay mode, the first and the second row are 
evaluated with and without the annihilation contributions, respectively. 

Table 1: Experimental data [2] and theoretical predictions for the branching ratios (in units 
of 10 -6 ). The numbers in columns Case I and Case II are our fitting results with gx and St 
constrained by varying the experimental data within la and 2a error bars, respectively. For 
each decay mode, the first (second) row is evaluated with(out) the annihilation contributions. 

Decay channel Experiment SM Case I Case II 

B-^r)K~ 2.2 ±0.3 2.21+J;! 2.29+°;^ 2.19+°;! 

1.80 ±0.82 2.19 ±0.27 2.17 ±0.34 
B°^ V K° <1.9 1.34+^ 1.30133 l-25±g^ 

1.02 ±0.65 1.42 ±0.22 1.34 ±0.30 
B-^rjK*- 19.5+H 6.27+j;^ 18.15+°;;* 17.44+°;! 

5.03 ± 1.82 17.88 ±0.49 17.45 ± 0.61 
B°^ V K*° 16.1 ± 1.0 6.87±1;1 16.96+°;°6 17.13±g;|| 

5.60 ±1.95 17.10 ±0.29 17.03 ± 0.65 



From Table (TJ we can see that the theoretical predictions within the SM for both B(B~ — > 
rjK~) and B(B° — > t]K ) agree with the experimental data within errors. However, both 
B(B~ — > 7]K*~) and B(B° — > rjK* ) are quite lower than the experimental data. We also note 
that our results are a little different with those in Refs. [EJ [HI [23], due to different choices for 
the input parameters, such as the moderate strength of Xa, Ab, and so on. 

As shown in Fig. [HI the four B(B°' + —> rjK^ ,+ ) are very sensitive to the presence of the 
7~Ce/f of Eq. ([H]). The two bands for the parameter spaces constrained by B~ —> r\K~ and 
B° — > r)K° decays are much overlapped, and the same situation is also found for the two bands 
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Figure 3: The contour plots for the parameters gx and 5t with the experimental data varying 
within la ((a) and (c)) and 2a ((b) and (d)) error bars, respectively. The red triangle, green 
triangle, blue circle, and yellow circle bands come from the decays B~ — > rjK~ , B° — > r)K°, 
B~ — ► rjK*~ , and B° — > r]K*°, respectively. Plots labels 'with(out) anni.' denote the results 
with(out) the annihilation contributions. 



constrained by B~ — > r]K*~ and B° — > rjK* decays. However, we note that B(B — > r/K) 
and BiB — > r]K*) have quite different dependence on these NP contributions. So, the allowed 
regions for the parameters gr and 5t are severely narrowed down when constraints from these 
four decay modes are combined. 

The final allowed regions for the parameters gx and 5t extracted from the four B 0,+ — > 
^/^(*)o,+ decays are shown in Fig. HI where the left (right) plot is the results with(out) the 
annihilation contributions. In addition, the dark and the gray regions in Fig. H] correspond to 
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without anni . 
1 "mPpiK**-': 



0.5 1 1.5 0.5 1 1.5 

g T [xlO~ 2 ] g T [xlO~ 2 ] 

Figure 4: The allowed regions for the parameters gx and 5t constrained by the four B°' + — > 
r]K^°' + decays. The dark and the gray regions correspond to the results obtained with the 
measured branching ratios varying within lcr and 2a error bars, respectively. The left (right) 
plot denotes the results with(out) the annihilation contributions. 




the results obtained with the measured branching ratios varying within la and 2a error bars, 
respectively. From now on, we denote these two possible regions by Case I and Case II. Com- 
paring the two plots in Fig. HI one can find that the uncertainties of annihilation contributions 
would loosen constraints on the NP parameters, especially on 5t- The numerical results for the 
parameters qt and 5t corresponding to the above two allowed regions are presented in Table EJ 
As shown in Table [T], with the parameters qt and 6t varying within these two allowed 
regions, the large B{B~ — > 7]K*~) and B(B° — > t]K*°) can be accounted for by the anomalous 
four-quark tensor operators without violating B(B~ — > r]K~). It is also interesting to note that, 
taking the 90% CL upper limit of B(B° — > t]K ) as an input, our fitting result for B(B° — * rjK°) 
is in good agreement with the very recent measurements 

B(B° -> 7]K°) = (1.8ig;g ± 0.1) x 10~ 6 BABAR , (27) 
B(B° ^7]K°) = (1.1 ±0.4 ±0.1) x 10~ 6 Belle [6], (28) 

which give the average value B(B° — > r/K ) = (1.3 ± 0.3) x 10~ 6 . 

3.2 Constraints on the NP parameters from B° —> (j)K* {) decay 

Now we discuss constraints on the NP parameters from B° — > <pK*° decay. With the NP weak 
phase St already extracted from the four B 0,+ — > r]K^ 0,+ decays, the parameter g' T could be 
severely constrained by the well measured B(B° —>■ <f)K*°) and fi- Our results are presented in 
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Table 2: The numerical results for the parameters qt and 6t, corresponding to the two allowed 
regions shown in Fig. UJ The allowed regions for the parameter g' T constrained by B° — > (j)K*° 
decay are also presented. For each case, the first (second) row denotes the results with(out) the 
annihilation contributions. 

Allowed region g r (xl0~ 3 ) g^(xl0 -3 ) 5y(rad) 
Case I dark 7.3j£!j 7.6iJ;g 0.77i°;?° 

7.7 ±1.6 8.9 ±0.5 0.70 ±0.13 
Case II gray 7.2±|;| 8.3^| 0.60±g;§| 

8.0 ±2.3 9.1 ±0.6 0.55 ±0.20 



Table 3: Experimental data [2] and theoretical predictions for the observables in B° — ► cpK* 
decay. The other captions are the same as in Table [TJ 

Observable Experiment SM Case I Case II 

B(xl0~ 6 ) 9.5 ±0.9 5.9±i;g 9.6l°| 9.91?;? 

5.7 ±0.6 10.0 ±0.3 10.3 ±0.7 

/l 0.49 ±0.04 0.761°;°? 0.50±°;°J 0.541°;°^ 

0.78 ±0.03 0.52 ±0.01 0.55 ±0.01 

0n(rad) 2.41±g:« 2.90i°;| 1.95^3 1.81±gS 

2.91 ±0.01 1.82 ±0.05 1.72 ±0.08 

± (rad) 2.52 ±0.17 2.90±°;^ L96 i°U2 i- 83 ^ 

2.91 ±0.01 1.82 ±0.05 1.73 ±0.08 



Table El Figs. El and El 

From Table El we can see that, with our choice for the input parameters, especially our quite 
optimistic choice for the annihilation contributions, the SM predictions for both B(B° — > <pK*°) 
and fi deviate from the experimental data, and possible NP scenarios beyond the SM may be 
needed to resolve the observed polarization anomaly. 

Fig. El shows the dependence of B(B° — > <pK*°) and fi on the parameter g' T , with the NP 
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Figure 5: The dependence of B(B° — > 4>K*°) and /l(-B° — ► ^if* ) on the parameter ^ with 
the NP weak phase St extracted from the four B°' + — > r]K^°' + decays. The upper and the 
lower plots denote the results with and without the annihilation contributions, respectively. In 
each plots, the solid blue (red) curves are the results with St given by Case I (II), and dashed 
curves due to the error bars of this parameter. The horizontal lines are the experimental data 
with the solid lines being the central values and the dashed ones the error bars (la and 2a). 



weak phase St extracted from the four B°' + — > decays. To illuminate the dependence 

more clearly, we have taken all the other input parameters to be their center values. As shown 
in Fig. [5j both B(B° — > (j)K*°) and fi are very sensitive to the parameter g' T . It is particularly 
interesting to note that the variation trends of these two observables relative to the parameter g' T 
are opposite to each other. Thus, with the allowed regions for the NP weak phase 5t extracted 
from the four B°' + —> i]K ( -* S)0 ' + decays, and the experimental data on these two observables, we 
could get constraints on the parameter g' T . The final allowed regions for the parameters g' T and 
5t are shown in Fig. [Ul with the corresponding numerical results given in Table [21 

Corresponding to the two allowed regions for the parameters g' T and 8t given in Table [2j 
both B(B° — > <f)K*°) and fi are in good agreement with the experimental data as shown in 
Table [HJ On the other hand, without the annihilation contributions, our predictions for the 
relative phases <p\\ and <fr± in the decay B° — > <pK*° are not consistent with the experimental 
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Figure 6: The allowed regions for the parameters g' T constrained by B(B° — > (j)K*°) and Jl 
with the NP weak phase 5t extracted from the four B°' + — > rjK^°' + decays. Other captions 
are the same as in Fig. UJ 



data. This mismatch is, however, reduced by including the annihilation contributions associated 
with strong phase, which may indicate the annihilation contributions to be complex. We 
also note that the annihilation contributions have been proved recently to be real and power 
suppressed [ID]. Moreover, in Ref. [JT], it is found that the NP strong phases are generally 
negligibly small compared to those of the SM contributions. Further discussion of possible 
resolution to the mismatch would be beyond the scope of this paper. 

3.3 Final allowed regions for NP parameters and predictions for 
B s — > decay 

Finally, using the relations 

9T = Ct8 + ^, (29) 
S'r = (l + ^r)Cr 1 + (i + ^)C T8 , (30) 

with N c = 3, we get the final allowed regions for the parameters CVi and CVs, which are 
presented in Table HI It is interesting to note that the color-octet operator Ots dominates the 
NP contributions. Actually, with T8 only, we obtain C T8 = (7.3±{ j) x 1(T 3 ((7.7±1.6) x 10~ 3 ) 
and C T8 = (9.1±i:g) x KT 3 ((10.7 ± 0.6) x 10" 3 ) from B -> r]K^ and B -> <j)K* decays, 
respectively, where (also for the following results) the numbers in the bracket are obtained 
without the annihilation contributions. Thus, we obtain a single color-octet 0^8 solution with 
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Table 4: Final results for the coefficients Cpi, Cps, and the weak phase 5t extracted from 
B — > 7]K^*> and B° — > (/)K*° decays. The other captions are the same as in Table 

Allowed region C T1 (xlO" 3 ) CV 8 (xl(T 3 ) 5 T (rad) 

Case I dark 1.7+];? 6.8±?;1 0.77±g;?g 

2.8 ±1.6 6.7 ±2.1 0.70 ±0.13 

Case II gray 2.6±|;g 6.4+|f 0.60j;{j;ijjj 

2.8 ±2.2 7.1 ±3.0 0.55 ±0.20 



C T8 = (8.5+°^) x 10~ 3 ((10.3 ±0.5) x 10~ 3 ) and 5 T = 44.1°™° (40.2° ± 7.2°). However, with 
the color-singlet operator Ot\ only, we could not get any solution. 

It is noted that the two tensor operators Oti and Ots could be expressed by the other two 
(pseudo-) scalar operators Os+p and 0' S+P through the Fierz transformation relations [IB] 

Os+p = t^Oti — -0 T 8, O's+p = TT^^ts — ~Oti- (31) 
12 6 12 6 

Such operators could be generated in many NP scenarios with scalar interactions. It would be 
useful to present the constraints on the coefficients of the two (S + P) <S> (S ± P) operators. To 
this end, we get 

C s+P = (-0.99+°;!) x 10 ~ 3 ((-0.89 ±0.38) x 10~ 3 ), (32) 
C' s+P = (0.29±g;H) x 10~ 3 ((0.09 ±0.32) x 10~ 3 ), (33) 

with the normalization factor ^§|Vts| and a new weak phase 5t = 44.1°+g 1 o (40.2° ± 7.2°), 
corresponding to the two tensor operators case. The single Ot% solution would correspond to 
C s+ p = -2C' S+P = (-1.4^) x 10~ 3 ((-1.7 ± 0.1) x 10~ 3 ), and the weak phase is the same 
as the former. It should be noted that the relation Cs+p = —2C' S+P is strictly required for the 
last correspondence. 

To further test the particular NP scenario with anomalous tensor operators Eq. dHJ), we also 
present our theoretical predictions for the branching ratio and the longitudinal polarization 
fraction of B s — > <jxj) decay, as summarized in Table [5j Within the SM, its decay amplitude is 
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Table 5: Theoretical predictions for B s — > (fxj> decay both within the SM and in the particular 
NP scenario Eq. (JHJ). The other captions are the same as in Table [TJ 



Observable 


Experiment 


SM 


Case I 


Case II 


i3(xl(T 6 ) 




20.61S 


29.4^ 


20.1±1;1 






17±2 


30 ±3 


29 ±3 


h 




n oq+0.04 


u - ' °-0.08 


n 79+ - 07 

u - '^-0.07 



0.83 ±0.02 0.66 ±0.03 0.67 ±0.03 



given by [15] 



n P,h , p,fr _ 1 o,h _ £ p,h , np.h _ ±np,h 
"3 ' ^4 2 3,EW r,"4,EW^- u 3 H,EW 



T/J4 2^ 4 ' EW 



(34) 



With respect to the relevant quantities in Eq. (J34J), one can get them directly from those for 
B° -> 0if* decay with some simple changes. This decay mode is of particular interest to test 
the proposed resolutions to the polarization anomalies observed in B — > (f)K* decays, since both 
of these decay modes are mediated by the same quark level subprocess b —* sss. In addition, 
since the two final-state mesons are identical in this decay mode, more observables in the time- 
dependent angular analysis will become zero [39J. So, this decay mode can be considered as an 
ideal probe for various NP scenarios proposed to resolve the polarization anomalies observed 
in B — > cf)K* decays. The earlier studies of this particular decay mode within the QCDF 
formalism have been carried out in Ref. [12], but without taking into account the annihilation 
contributions. 

From the numerical results presented in Tables [3] and [51 we can see that both the branch- 
ing ratio and the longitudinal polarization fraction of B s — > cfxfi decay is larger than those 
of the decay B° — > <pK*°, which is due to the relative factor of two in the cfxp amplitude 
Eq. (134j) . as well as the additional contribution from the annihilation coefficient /3^ h . In 
particular, due to an accidental cancelation, the annihilation coefficient (3 P '° (contributing to 
both H 00 (B° —> 4>K*°) and H 00 (B S — > (jxp)) is quite smaller than [T3]. Interestingly, 
recent calculations made in the perturbative QCD approach also predict a large branching 
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ratio B(B S —> (pep) = (44. l^: 3 ^ .4^.0) x 10 ~ 6 witn tne longitudinal polarization fraction 
f L (B s -> 00) = (68.0j£gj£££{j) x !0" 2 P]. Their result for f L (B s -> 00) is smaller than our 
predictions within the SM. The predicted results presented in Table [5] could be tested at the 
Fermilab Tevatron and the LHC-b experiments. 

4 Conclusions 

In summary, motivated by the observed discrepancies between the experimental data and the 
SM predictions for the branching ratios of B — > r/K* decays and the polarization fractions in 
B — > cf)K* decays, we have studied a particular NP scenario with anomalous tensor operators 
Oti and Ots, which has been proposed to resolve the polarization anomalies in the literature [T71 

EEUDl. 

After extensive numerical analysis, we have found that the above observed discrepancies 
could be resolved simultaneously and constraints on the NP parameters have been obtained. 
With both the experimental data and the theoretical input parameters varying within la error 
bars, we have found the following two solutions: (I) both the two tensor operators contribute 
with the parameter spaces presented in the upper two rows of Table HI (II) only Ots contributes 
with C T 8 = (8.5+o;g) x 10" 3 ((10.3 ± 0.5) x 10~ 3 ) and S T the same as solution (I). The above 
two solutions correspond to the scenario with new operators Os+p and 0' S+P added to the 
SM, with the parameter spaces: (I) C s+P = (-0.99+°:|) x 10~ 3 ((-0.89 ± 0.38) x 10" 3 ), 
C' s+P = (0.29l{j;ll) x 10- 3 ((0.09 ± 0.32) x 10~ 3 ) and (II) C s +p = -2C' S+P = (— 1.4±g:|) x 
10~ 3 ((— 1.7 ± 0.1) x 10~ 3 ). Their weak phase 5s+p is the same as St- 

Against our early expectations, we have found that the solution with two tensor operators 
are dominated by the color-octet operator Ots- It would be interesting to investigate whether 
the available NP models could give such effective interactions at the m;, scale. 

To further test the particular NP scenario with two anomalous tensor operators Oti,ts 01 
(pseudo-) scalar operators 0^ +p , we have also presented our predictions for the observables 
in B s 00 decay, which could be tested more thoroughly at the Fermilab Tevatron and the 
LHC-b experiments in the near future. 

It should be noted that the strong constraints in Table H] obtained without annihilation 
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contributions may be too optimistic. As has been shown in the table, the uncertainties due 
to poorly known annihilation contributions could dilute the requirement of NP very much. 
Generally, this caveat could be applied to probe possible NP scenarios in exclusive non-leptonic 
B decays. Further theoretical progress is strongly demanded. 
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Appendix A: Decay amplitudes in the SM with QCDF 



The amplitudes for B —> rjK^*' are recapitulated from Ref. 



Yl V pb V P*s\ A K-r, q S pu a 2 + 2o| + -a p 3 

p=u,c 

+V2A 



EW 



c n | p | p p P i aP i aP 

Opu P2 + « 3 + « 4 _ 2 a 3' EW ~ 2" 4 ' EW ^ 3 ft,EW 



+A 



5pu («1 + P2) + «4 + «4,EW + 03 + 03,EW 



(35) 



B — *K r] 



X VpbV* a <j Aj{ t 
p=u,c 

+V2Af< 



8pu &2 + 2«3 + -a 



3,EW 



1 



as + a A — —a 



2"3,ew 2 aA ' EW 



+ 03 2 03,~E\ 



a 4 ~~ 2" a 4,EW 



+ 03 ~ -ff 



3 2 M 3,EW 



(36) 



\/2A SM 



p=u,c 

+V2A 



' 'hi 



b pu oi2 + 2af + -a 



p 

3,EW 



r a 1 p 1 p ^ P ^~ P 1 OP 1 /3P 

<VP2 + a 3 + « 4 - -a 3i Ew ~~ 2° 4 ' EW ^ 3 ^,ew 



+A 



tfpu («1 + 02) + «4 + «4 EW + 03 + P3 EW 



(37) 
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B — >K r] 



E VpbVps \ AR*o Vq 

p=u,c 



-a 



p 

3,EW 



oR-\-a p A a 



3,EW 



-a 



V 



4,EW 



+ A? - M 



-A 



j a 4,EW + 03 2^3,EW 



(38) 



where the explicit expressions for the coefficients ccf = af (MiM 2 ) and /3f = flf (MiM 2 ) could 
also be found in Ref. [23]. We recall that the terms contain one-loop vertex, penguin and 
hard spectator contributions, whereas the f3f terms are due to the weak annihilation, and the 
transition form factors are encoded in the factors Am 1 m 2 - 

The decay amplitude of B — > <pK mode can be read off from Refs. [T5l [T6] 



^S^TT - Ah K*°4> E V p* V ps 



p,h . p,h 



p,h 



a 3 ' + a I a 3,EW n a 4,EW 



+ n o h,ew 



(39) 



p=u,c 



where the superscript '/i' denotes the helicity of two final-state vector mesons, with h = 0, +, — 
corresponding to two outgoing longitudinal, positive, and negative helicity vector mesons, re- 
spectively 

In the helicity basis, the decay amplitude Eq. (1391) can be further decomposed into the 
following three helicity amplitudes 



ttSM 
^00 



H 



SM 



with us 



^ E W 

p=u,c 

4-, E w 



p,0 , p,0 
«3 + «4 



-a 



p,0 
3,EW 



a 



p,0 
4,EW 



p,0 



,n,Ew 



„p.± 



2 



tt 3 ' + a 4 ' a 3,EW o a 4,EW 



1 

— ( 

2 



"r ^3 2^ 3 ' EW 



p=u,c 



A h R ,o, = ^ (K*°|(^)^|fi°)(0|(^)y|O) 



(40) 
(41) 

(42) 



4° 



A ± 



iG F m Bd U 
y/2 2m K *o 
iGp 



m Bd 
rriK*° 



B-*K* 



1 + 



m Bd 



-)A\ 
m Bd 

rriK*° 
m Bd 



[m. 



B->K* 



(43) 
(44) 
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Appendix B: Theoretical input parameters 



Bl. Wilson coefficients and CKM matrix elements 

The Wilson coefficients Ci(p) have been evaluated reliably to next-to- leading logarithmic or- 
der pHIHj. Their numerical results in the naive dimensional regularization scheme at the scale 
H = m b (fi h = ^/A h m b ) are given by 

Ci = 1.077 (1.178), C 2 = -0.174 (-0.355), C 3 = 0.013 (0.027), 

C 4 = -0.034 (-0.060), C 5 = 0.008 (0.011), C 6 = -0.039 (-0.081), 

C 7 /a e . m . = -0.013 (-0.034), C 8 /a e . m . = 0.047 (0.099), C 9 /a e . m . = -1.208 (-1.338), 

Cio/cw = 0.229 (0.426), C 7l = -0.297 (-0.360), C 8g = -0.143 (-0.168). (45) 

The values at the scale ph, with m b = 4.79 GeV and A h = 500 MeV, should be used in the 
calculation of hard-spectator and weak annihilation contributions. 

For the CKM matrix elements, we adopt the Wolfenstein parameterization [45] and choose 
the four parameters A, A, p, and rj as 

A = 0.809 ± 0.014, A = 0.2272 ±0.0010, 7J = 0.197+°;°^, fj = 0.339^;°^, (46) 

with p — p (1 — 4r) and 77 = 77 (1 — 4r). 

B2. Quark masses and lifetimes 

As for the quark mass, there are two different classes appearing in our calculation. One type is 
the pole quark mass appearing in the evaluation of penguin loop corrections, and denoted by 
m q . In this paper, we take 

m u = md = m s = 0, m c = 1.64 GeV, m b = 4.79 GeV. (47) 

The other one is the current quark mass which appears in the factor r^f through the equation 
of motion for quarks. This type of quark mass is scale dependent and denoted by fn q . Here we 
take [1J 

m s (p)/m q (p) = 25 ~ 30, m s (2GeV) = (95±25)MeV, m b (m b ) = 4.20 GeV, (48) 
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where m 9 (/x) = (fn u + m d )(/i)/2, and the difference between u and d quark is not distinguished. 

As for the lifetimes of B mesons, we take pQ tb u = 1.638 ps, rs d = 1.530 ps, and tb s = 
1.466 ps as our default input values. 

B3. The decay constants and form factors 

In this paper, we take the decay constants 

f B = (216 ± 22) MeV [47], f Bs = (259 ± 32) MeV [47], f w = (130.7 ± 0.4) MeV [I], 
fx = (159.8 ± 1.5) MeV P fx* = (217 ± 5) MeV [48], f^ = (231 ± 4) MeV [48], 
f K ,{l GeV) = (185 ± 10) MeV [IE], f£(l GeV) = (200 ± 10) MeV [IE], (49) 

and the form factors [45] 

Fq^ v (0) = 0.275 ± 0.036, F^ K (0) = 0.331 ± 0.041, A^ K * (0) = 0.374 ± 0.033, 

(0) = 0.411 ± 0.033, Af^ K * (0) = 0.292 ± 0.028, A^ K * (0) = 0.259 ± 0.027, 

V Bs ^(0) = 0.434 ± 0.035, Af s ~^(0) = 0.311 ± 0.030, Af fl ^(0) = 0.234 ± 0.028, 
if (0) = 0.333 ± 0.028, T^ K * (0) = 0.333 ± 0.028, if fl ^(0) = 0.349 ± 0.033, 
T 2 Bs ^ (O) = 0.349 ± 0.033. (50) 

For the parameters related to the 77-77' mixing, we choose [28] 

/, = (1.07 ± 0.02) U, f s = (1.34 ± 0.06) U, = 39.3° ± 1.0° . (51) 

The other parameters relevant to the meson 77 can be obtained from the above three ones [5]. 

B4. The LCDAs of mesons. 

For the LCDAs of mesons, we use their asymptotic forms [3UI EE] 

Ms) = ®l±( x ) = 9 { l )V (x) = 6x(l - x) , %{x) = 1, 

<&,(*) = 3 (2x - 1) , <# )V (x) = ^ [1 + (2x - l) 2 ] • (52) 

As for the B meson wave functions, we need only consider the first inverse moment of the 
leading LCDA $f (£) defined by [22] 

ff*f«)^, (53) 
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where A# = (460 ± 110) MeV [IH] is the hadronic parameter introduced to parameterize this 
integral. 
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